Effect of Delocalized Vortex Core States on the Specific Heat of Nb 
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The magnetic field (B) dependence of the electronic specific heat for a simple BCS type-II super- 
conductor has been determined from measurements on pure niobium (Nb). Contrary to expectations, 
the electronic specific heat coefficient j(T, B) is observed to be a sublinear function of B at fields 
above the lower critical field H c i- This behavior is attributed to the derealization of quasiparticles 
bound to the vortex cores. The results underscore the ambiguity of interpretation that arises in 
specific heat studies of this kind on newly discovered type-II superconductors, and also emphasize 
the need to do such measurements under field-cooled conditions. 

PACS numbers: 



Specific heat measurements on type-II superconduc- 
tors in the vortex state are sensitive to low-energy quasi- 
particle excitations. In fully-gapped superconductors the 
low-lying excitations are traditionally thought of as be- 
ing confined to the vortex cores 0. In this case the 
coefficient of the linear term in the specific heat c as a 
function of temperature, i.e. the 'Sommerfeld coefficient' 
7, is proportional to the density of vortices, and hence a 
linear function of the internal magnetic field B 2] . The 
situation is different in superconductors with a highly 
anisotropic energy gap or gap nodes. In these cases there 
is a major contribution to the electronic specific heat at 
low temperatures from delocalized quasiparticles. The 
flow of superfluid around a vortex lowers the energy of 
quasiparticles delocalized near the gap minima, resulting 
in a finite density of states (DOS) at the Fermi energy 
and a corresponding nonlinear field dependence for 7(B) 
3]. For an anisotropic gap, such behavior is expected 
only if the energy shift exceeds the gap minimum. Thus, 
7(B) tells us something about the gap structure. 

What has not been considered in the interpretation of 
many such studies, is that even for the case of a simple 
BCS superconductor characterized by an isotropic energy 
gap, the quasiparticles are not strictly confined to single 
vortices. Quasiparticle states bound to the vortex cores 
become delocalized when their spatially extended wave 
functions overlap those from nearest-neighbor vortices. 
The overlap is enhanced at higher magnetic fields where 
the vortices are closer together, but there is also some 
overlap as soon as there is more than one vortex in the 
sample. Likewise, the degree of derealization increases 
with increasing temperature due to thermal population 
of higher-energy core states, characterized by wave func- 
tions extending further out from the core center. This is 
the so-called 'Kramer- Pesch (KP) effect' 0. 

Recently, the field dependence of the vortex electronic 
excitations has been confirmed by thermal conductivity 
[f| and muon spin rotation (/iSR) Q experiments on 
VaSi. By detecting the change in the vortex core size, 
/iSR measurements on have confirmed that the de- 
localized quasiparticles detected by thermal conductivity 



originate from the vortex cores. These delocalized core 
states should have a profound influence on the field de- 
pendence of the specific heat, which probes the electronic 
DOS at the Fermi level, both inside and outside of the 
vortex cores. Recent calculations Q in the framework of 
the quasi-classical Eilenberger theory show that for an 
isotropic energy gap, 'y(B) ccB behavior persists only up 
to a crossover field B* , above which the overlap of core 
states results in a sublinear dependence of 7(B) on B. 
Strictly speaking B* is temperature dependent, reduced 
by thermal population of the more energetic and spatially 
extended core states. 

To date there has been no clear experimental verifica- 
tion of a crossover field B*(T) in a fully-gapped type-II 
superconductor. While there have been numerous spe- 
cific heat studies on BCS type-II superconductors, rarely 
is 7(B) observed to be a linear function of B. One rea- 
son is there are very few "simple" BCS type-II super- 
conductors. For example, specific heat measurements 
on the widely studied superconductor 2iJ-NbSe2 show 
that l(B) is a nonlinear function of B even at low fields 
IE E3 ■ There is now good evidence that 2iJ-NbSe2 is 
a multiband superconductor, with a different size energy 
gap on two different sheets of the Fermi surface pj [Tl| • 
The nonlinearity of i(B) can be attributed to two-band 
superconductivity Il2| . giving rise to a large vortex-core 
size at low fields [3, LL3] • 

Even the cubic A-15 compound VaSi is not so simple. 
For a magnetic field applied along the (001) direction, 
there is a substantial range of fields over which the hexag- 
onal vortex lattice either gradually distorts into a square 
lattice 0, ^| or coexists with a square lattice ^| . This 
vortex lattice transition is accompanied by derealization 
of quasiparticle core states PJ , and occurs at very low 
fields. Consequently, the specific heat is a linear function 
of B only over a very narrow field range far below the 
upper critical field H C 2 [T^ |. 

In view of these complications, here we investigate the 
field dependence of j(B) in the elemental type-II super- 
conductor niobium (Nb). A (OOl)-oriented, 99.99 % pu- 
rity Nb single crystal was obtained from Goodfellow Cor- 
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FIG. 1: Temperature dependence of the specific heat at sev- 
eral different applied magnetic fields, plotted as c(T, H)/T vs 
T 2 . The measurements were taken under FC conditions. The 
solid line through the normal state data at H — 5 kOe is a 
fit described in the main text. Inset: Extrapolation of the 
low-temperature data to K. 



poration. The sample was in the shape of a short cylin- 
der, ~ 1.8 mm long with a diameter of ~2 mm. Magne- 
tization measurements were made with a SQUID mag- 
netometer and specific heat measurements were made 
using a Quantum Design Physical Properties Measure- 
ment System (PPMS) that utilizes a thermal relaxation 
calorimeter. The calorimeter is composed of a sapphire 
sample holder and the addenda heat capacity is essen- 
tially field-independent for the fields employed in this 
work. The external field H was directed along [111], 
which under field-cooled (FC) conditions produces a well- 
ordered hexagonal vortex lattice in Nb for all fields H c i < 
H <H c2 [13. FC and zero-field cooled (ZFC) swept-field 
measurements of the sample magnetization M(T, H) and 
the specific heat c(T, H) were found to be greatly ham- 
pered by flux jumps. Consequently, the FC measure- 
ments reported here were done by field cooling the sam- 
ple from above T c to the desired temperature for each 
value of H . On the other hand, the ZFC measurements 
were done by first cooling the sample in zero field and 
then increasing H . From ZFC measurements of M(T) at 
5 Oe, the Nb crystal was determined to have a sharp su- 
perconducting transition temperature of T c = 9.2 K. The 
zero-temperature extrapolated value of the upper critical 
field was determined to be -ff C 2(0)«4.5 kOe, from which 
the superconducting coherence length is calculated to be 
£ = 270 A. Just above T c the resistivity of the sample 
saturates with a value of 0.32 /iQ-cm, from which we cal- 



culate the electron mean free path to be I = 2300 A. Thus, 
our sample is well within the clean limit. 

Measurements of the sample magnetization M (H) 
were done to determine the internal magnetic field 

B = H - AttM (1 - N) , (1) 

where N = 0.379 is the demagnetization factor of our 
Nb crystal. This value of N was determined by ZFC 
measurements of M(H) at several temperatures in the 
Meissner phase, assuming full shielding (B — 0). Under 
FC conditions the quantity AttM was found to be very 
small (with a maximum value of 31 G at H = 750 Oe and 
T = 2.25 K), indicating that very little flux is expelled. 
Since BmH, it is safe to interpret our measurements of 
7 (T, H) as 7 (T, B). We note that this equivalence is often 
assumed in specific heat studies without justification. 

The normal-state specific heat of Nb is described by 
the relation 

c„(T) = 7n T + /3T 3 , (2) 

where j n is the Sommerfcld constant in the normal state, 
and 7 n T and /3T 3 are the electronic and phonon contri- 
butions, respectively. FC measurements of the tempera- 
ture dependence of the specific heat of Nb are plotted 
in Fig. 1 as c(T,H)/T against T 2 . The solid line in 
the main panel is a fit of the H = 5 kOe data below 
10 K to Eq. ©, yielding 7n =7.60±0.02 mJ/mol K 2 and 
/3 = 0.14±0.01 mJ/mol K 4 . 

At T<^T C the specific heat is described by 

c(T,H) = 7(r,i/)T + /3r 3 + Ces(T,0) (3) 
= j(T,H)T + c(T,0), (4) 
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FIG. 3: (a) Temperature dependence of [c(T, H) -c(T, 0)]/T 
at several fields measured under FC conditions, (b) Field 
dependence of [c(T, H) — c(T, 0)]/T at several temperatures 
below T c . The open and solid symbols correspond to ZFC 
and FC measurements, respectively. The downturn at high 
fields for each data set indicates H C 2(T). For visual clarity 
the data at T= 3.0, 4.5 and 5.0 K are vertically offset. 



where c cs (T, 0) is the electronic specific heat in the su- 
perconducting medium surrounding the vortex cores. 
To determine the Sommcrfcld coefficient 7 (if) = 
liiriT-»o c(T, H)/T, the low-T data shown in the inset of 
Fig. 1 were extrapolated to T — * K assuming 



c(T,H) 
T 



7(T, H) + (3T 2 



-bT a /T 



T 



(5) 



The last term is the limiting low-temperature BCS ex- 
pression for c es (T, 0) divided by T. At If = kOe, the fit 
to Eq. © yields c es (T, 0)/j n T c = 10.19 exp(-1.64T c /T). 
As shown in Fig. 2, 7(B) is a linear function of B only 
up to H^0.22H c2 . 

The temperature and magnetic field dependences of 
the difference [c(T, H)-c(T,0)]/T are plotted in Fig. 3. 
At low temperatures this quantity reflects j(T, H), as in- 
ferred from Eq. In Fig. 3(a) we see that below 2.5 K, 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 

H (kOe) 



FIG. 4: Magnetization and specific heat data at T — 2.3 K. 
(a) FC (solid circles) and ZFC (open circles) measurements 
of dy/dH, and (b) ZFC measurements of dM/dH. The ZFC 
measurements were done by first cooling the sample in zero 
field to T = 2.3 K and then increasing H. 



[c(T,H)-c{T,0)]/T decreases with decreasing T. We 
attribute this behavior to shrinking of the vortex cores 
(i.e the KP effect), which reduces the contribution of the 
zero-energy DOS per vortex to the specific heat. The ef- 
fect is more pronounced at higher H, due to the increased 
density of vortices. Above 2.5 K, [c(T, H)-c(T, 0)]/T de- 
creases with increasing T due to the subtraction of the 
H = specific heat jump at T c , and hence does not re- 
flect the temperature dependence of 7. From the field 
dependence of [c(T,H) — c(T,0)]/T measured under FC 
conditions (see Fig. 3(b)), we conclude that j(B)(xB up 
to a temperature-dependent crossover field. 

To clearly see the field range over which 7 is a linear 
function of B, in Fig. 4(a) we show the derivative d^/dH 
of the 7« [c(T, H)-c{T, 0)]/T data at T = 2.3 K. For both 
FC and ZFC measurements, d'y/dH is constant only at 
low fields. In this same low-field range, j(T, H) is deter- 
mined to be zero from the ZFC measurements. Thus, the 
nonzero value of d^ /dH in the FC data at low fields must 
result from trapped vortices below the lower critical field 
H&. To see that this is indeed the case, we determined 
the field of first vortex entry H cn (T) from ZFC measure- 
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ments of M(H) for increasing H . As shown in Fig. 4(b), 
the increase of dM/dH at H en agrees with the decrease 
(increase) of drf/dH in the FC (ZFC) measurements of 
the specific heat. This was also determined to be the 
case at the other temperatures. We note that H c i(T) is 
somewhat lower than H e n(T), if one assumes there is a 
Bean-Livingston barrier that both impedes the en- 
try of flux when the field is increased above H C \(T), and 
traps flux below H C \(T) in the FC measurements. 

The observation of the simple relation 7(B) cx B 
only in FC measurements below H ea (T), indicates that 
the trapped vortices form a highly-disordered lattice, 
in which they behave as isolated vortices. The disor- 
der apparently disrupts the connection between the lo- 
cal DOS of nearest-neighbor vortices that is found in 
a well-ordered lattice |l9j. At fields immediately above 
H en (T), 7(-E>) exhibits a sublinear dependence on B. As 
explained in Ref. the sublinear dependence of j(B) 
on B is due to the shrinkage of the vortex cores that oc- 
curs as a result of the derealization of the higher-energy 
core states. While this may seem at odds with electronic 
thermal conductivity measurements of Nb that indicate 
only a small degree of quasiparticle derealization just 
above H c i |2Cj , the details of the experimental method 
must be considered. In Ref. j2£j the thermal conductiv- 
ity k(H) was measured in both monotonically increasing 
and decreasing field. In Fig. 3 the increasing field ZFC 
measurements of c e (T, H) exhibit an upward curvature 
above H en . This indicates that the flux has some diffi- 
culty entering the sample even above H on . We note that 
the measurements of Ref. [2(j are often cited in the liter- 
ature as the most archetypal example of n{H) for a clean 
conventional type-II superconductor. However, here we 
see that the intervortex transfer of quasiparticles at low 
fields in Nb is greatly enhanced for a highly-ordered vor- 
tex lattice generated under FC conditions. 



The crossover field for an isotropic-gapped supercon- 
ductor was determined in Ref. to be B*(Q.l T c ) w 
0.337J c2 . Here we find for Nb that B*(0.1 T c ) < 
-ff cn (0.1 T c ) «0.21i7 C 2. The ratio of the minimum to max- 
imum superconducting energy gaps on different Fermi 
sheets in Nb is estimated to be ~0.8 21]. According to 
Ref. [3| , gap anisotropy of this size reduces the crossover 
field to £*(0.1 T c ) fa0.25H c2 . Thus, the small difference 
in gap values may be sufficient to explain the complete 
absence of a crossover field in the vortex state of Nb. 

In summary, the T-linear coefficient of the electronic 
specific heat j(T, B) has been determined for pure Nb 
in the vortex state. Contrary to popular belief, 7(T, B) 
is not a linear function of B at fields immediately above 
H c \. Our results support theoretical calculations show- 
ing that even in a simple BCS type-II superconductor, 
7(T, B) is a sublinear function of B over a wide region of 
the vortex phase. This calls into question specific heat 
studies on type-II superconductors that have reported 
l{B) cx B behavior persisting up to fields close to H c2 . 
We suspect inaccuracies in the T — > K extrapolation 
procedure are partially responsible, but also stress the 
importance of small field increments in specific heat mea- 
surements of this kind. Lastly, the effect of delocalized 
quasiparticle core states on the field dependence of the 
specific heat should be considered when attempting to 
identify the pairing symmetry of a newly discovered su- 
perconductor by this method. In other words, the ab- 
sence of j(B) oc B behavior is not necessarily an indica- 
tion of unconventional superconductivity. 
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